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1 . INTRODUCTION 


Much  progress  has  been  made  in  recent  years  in  developing  viable 

deterministic  models  of  crack  growth  in  gun  tubes.  Such  models  relate 

the  eunount  of  growth  during  a given  number  of  cycles  to  various 

material  properties  and  design  parameters.  Using  such  models  one  can 

then  predict  the  number  of  cycles  required  for  the  tube  to  fail, 

provided  one  knows  all  the  relevant  material  and  design  properties. 

In  practice,  such  properties  vary  from  gun  to  gun  and  may  even 

vary  from  round  to  round  throughout  the  life  of  a single  gun.  Gun 

tube  fatigue,  therefore,  cannot  be  treated  as  a purely  deterministic 

process.  It  is  necessary  to  allow  a place  for  random  elements  which 

must  be  analyzed  statistically  rather  than  deterministically.  The 

Army  has  recognized  this  fact  by  stating  gun  tube  fatigue  requirements 

in  statistical  terms. 

1 2 

Racicot  and  Proschan  and  Sethuraman  have  devised  probabilistic 
models  of  gun  tube  fatigue  by  adding  stochastic  elements  into  an 
existing  deterministic  model  of  fatigue.  It  is  the  aim  of  this  paper 
to  examine  this  probabilistic  approach  in  somewhat  more  detail  and  to 
examine  somewhat  more  closely  the  effects  of  various  assumptions  and 

^Racicot,  R.  L.,  "A  Probabilistic  Model  of  Gun  Tube  Fatigue," 
Watervliet  Arsenal  Technical  Report  No.  ARLCB-TR-77029  Watervliet,  NY 
May  1977.  All  further  references  to  Racicot  are  to  this  report. 

2 

Proschan,  F.,  and  Sethuraman,  J.,  "A  Probabilistic  Model  for  Initial 
Crack  Size  and  Fatigue  Life  of  Gun  Barrels,"  Florida  State  University 
Technical  Report  No.  M394,  The  Florida  State  University,  Tallahassee, 
Florida,  December  1976. 
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methods  of  these  authors  on  their  conclusions.  Racicot  has  already 

stated  the  benefits  to  be  derived  from  a closer  study  of  the 

probabilistic  aspect  of  gun  fatigue;  we  shall  not  repeat  them  here. 

2.  A FRACTURE  MECHANICS  MODEL  OF  CRACK  GROWTH  CONSIDERED  STOCHASTICALLY 

The  deterministic  model  mentioned  above  was  developed  by  Throop 
3 

and  others  at  Watervliet  Arsenal  based  on  the  Paris  expression  for 
rate  of  crack  growth  and  experimental  results.  The  rate  of  crack 


growth  is  approximated  by: 


^ - (AK) 
dN  ~ M 


where:  b = crack  depth 

AK  = range  of  variation  of  stress  intensity  factor  K 
over  one  cycle 
N = number  of  cycles 

m = a parameter  dependent  on  material  properties  and 
stress  intensity.  Usually  m is  between  2 and  4. 
M = a parameter  dependent  on  material  properties. 

In  this  model,  AK  and  M are  given  by: 

AK  = aS/ifb 


= EKjg(Jy/C 


where : 


S = maximum  hoop  stress  at  the  bore 

a = a parameter  depending  on  crack  shape  and  the  residual 
stresses  introduced  by  the  autofrettage  process 

throop,  J.  F.,  and  Miller,  G.  A.,  "Optimum  Fatigue  Crack  Resistance," 
ASTM  Special  Technical  Publication  467,  Philadelphia,  PA,  1970, 
pp.  154-168. 


E = Young's  modulus 

Kj^  = fracture  toughness  for  a crack  in  a tangential 


stress  field 
cfy  = yield  strength 

C = a psirameter  varying  with  m to  maintain  dimensional 
homogeneity  and  possibly  depending  on  material 
properties. 

Combining  (1),  (2)  and  (3)  gives: 

db  ~k+l 

ST  ■ “ (4) 

where:  k = ^ - 1 

G - (5) 

EOyKlc 

Racicot  assumes  that  G and  k are  essentially  constant  during  tube 
life,  so  that  he  can  solve  the  differential  equation  (4)  by  simple 


integration  to  give  : 


1 -k  -k 

N = icbo’^-b^ 


where  bg  is  an  initial  crack  depth,  assumed  to  be  present  after  a few 
rounds  of  firing. 

A different  line  of  reasoning,  not  assvuning  the  constancy  of  G, 
can  also  be  used.  (This  line  of  reasoning  is  used  by  Cramer  in  a 
similar  problem.)  Assume  that  k is  essentially  constant  over  tube  life. 
Equation  (4)  may  be  rewritten: 

*Racicot  actually  used  N-Nj^  where  we  have  used  N.  is  the  number  of 

cycles  required  to  crack  initiation.  In  practice  is  small  enough 

to  be  ignored. 

^Cramer,  H.,  Mathematical  Methods  of  Statistics,  Princeton  University 
Press,  Princeton,  NJ,  1946,  p.  219. 


(7) 


- "j-1 


k+1 

= 


where  bj  = the  crack  length  after  the  jth  cycle. 

Gj  = the  value  of  G inunediately  preceding  the  jth  cycle. 

From  equation  (7),  we  have:  b--b-  ■, 

G.  = ^ 

J , k+1  (8) 


Summing  over  j gives: 


N N b.-b.  1 

3-1^  j=i 


(9) 


If  bj^j-bj  tends  to  be  relatively  small,  then  the  sum  on  the  right- 
hand  side  of  (9)  can  be  approximated  by: 

J -T-T  = r C^o  -^N  5 

b„  uk+1  k 


’o  b'' 

And  so  at  least  approximately: 


N 


1 


-k  u -k. 


j=l 


(10) 


Assuming  that  G is  constant  over  tube  life  means  that  Gj  = G for  all 

j and  hence:  N 

iGj  = 

j=l^ 


NG 


(11) 


Combining  (10)  and  (11)  then  gives  equation  (6)  for  constant  G. 

From  a probabilistic  viewpoint,  equation  (10)  is  very  nice. 
First,  it  shows  that  crack  growth  depends  probabilistically  on  four 
things: 


a.  The  distribution  of  the  initial  crack  depth  bg. 

b.  The  distribution  of  the  rate  exponent  k. 

c.  The  distribution  of  Gj,  which  is  the  value  of  G immediately 
following  the  forming  of  the  initial  crack  and  preceding  the  first  cycle. 

d.  The  manner  in  which  G probabilistically  varies  from  cycle  to  cycle. 


4 


Second,  equation  (10)  involves  a sum  of  random  variables  which 


% 


suggests  that  if  the  Gj's  are  sufficiently  well-behaved,  some 
classicial  results  such  as  the  central  limit  theorems  may  be  exploited. 

Also  note  that  the  argviment  used  to  derive  (10)  does  not  depend 
in  any  way  on  the  exact  material  properties  or  design  variables  that 
go  into  G,  nor  does  it  require  any  independence  (mathematical, 
physical  or  statistical)  of  G and  b. 

The  authors  cited  above  all  assume  that  G does  not  vary  from 
cycle  to  cycle.  (Racicot  did  allow  S to  vary  for  some  of  his 
calculations.)  We  shall  consider  the  effect  of  this  assumption  next. 


3.  EFFECTS  OF  CYCLE-TO-CYCLE  VARIATIONS  IN  THE  MODEL  PARAMETERS  ] 

There  is  no  reason  to  suppose  that  the  Gj's  are  actually  constant 

from  cycle  to  cycle.  However,  if  the  probability  distribution  of 
N 

^ G.  is  close  to  that  of  NGi,  in  some  suitable  sense,  then  NG,  can 
j=l  ^ N 

effectively  replace  ^ G.  in  equation  (10). 

j=l  ^ 

Since  fatigue  failures  typically  occur  only  after  thousands  of 
cycles,  it  is  quite  legitimate  to  assume  that  in  equation  (10),  N is 
large.  This  allows  us  to  make  good  use  of  the  central  limit  theorem 
and  the  Laws  of  Large  Numbers. 

We  will  examine  the  probabilistic  behavior  of  bj^  for  large  N 
under  some  assumption  as  to  how  Gj  varies  from  cycle  to  cycle.  For 
some  of  these  assumptions,  we  will  show  that  the  asymptotic  behavior 
(that  is,  the  behavior  as  N approaches  infinity)  is  not  much  different 
from  the  behavior  when  the  Gj's  are  assumed  not  to  vary  from  cycle  to 
cycle.  In  other  cases,  the  asymptotic  behavior  can  be  quite  different. 
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The  various  sets  of  assumptions  listed  below  are  not  meant  to  be 
exhaustive.  They  were  selected  simply  because  they  are  the  sort  of 
assumptions  that  spring  most  readily  to  a statistician's  mind.  They 
are  also  the  sort  of  assumptions  that  are  used,  often  without 
conscious  realization,  by  non-statisticians  who  employ  "cook-book" 
formulas  from  the  standard  texts.  In  the  discussion  below,  we  will 
employ  "error  terms"  denoted  by  ei,e2-..  We  assume  that  these  are 
independent  and  identically  distributed  with  common  mean  and 
common  standard  deviation  a^.,  where  0 < < <».  We  also  let  and 

Oq  be  the  mean  and  standard  deviation,  respectively  of  G^.  Let 

+ - + Gj^.  Then,  the  problem  is  to  determine  the  asymptotic 
behavior  of  - NGj^. 

Case  I:  Gj  = Cj 

This  assumption  implies  that  Gj,G2...  are  statistically  independent. 
However,  since  GpG2...  are  functions  of  material  and  design  parameters, 
which  should  not  change  that  radically  from  cycle  to  cycle,  the 
independence  of  the  G's  would  seem  unlikely.  We  shall  not  consider 
this  case  further. 

Case  II:  Gj  = Cj  + ^ for  j > 1 

This  is  a simple  linear  regression  model.  In  this  case, 

4 

Sn  = NGj  + 6]^  +....+  En-I*  Laws  of  Large  Numbers  then  imply 

that  (Sj^  - NG2^)/(N-1)  approaches  y^  in  probability.  In  particular  if 
Pg.  = 0,  then  distribution  of  may  be  approximated  by  that  of  NGj^. 

^Cramer,  H. , Mathematical  Methods  of  Statistics.  Princeton  University 
Press,  Princeton,  NJ,  1946,  p.  497. 
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How  good  the  approximation  will  be  depends  on  In  fact  it 

^ follows  from  the  Central  Limit  Theorem  that  for  large  N, 

1/2 

Sfj  = NGi  + (N-l)y^  + (N-1)  O^Un 
where  Uj^  has  approximately  a standard  normal  distribution.  When 
Pg.  = 0 and  Og.  is  sufficiently  small,  the  deviation  Sj^  - NGj  will  tend 
to  be  close  to  zero  for  large,  but  not  too  large  N. 

Case  III:  Gj  = for  ra  > 1. 

In  this  case,  Sj^  - NG]^  = (ei  + ...  + ej^_2-N)G2.  The  Laws  of 
Large  Numbers  imply  that  (S|^  - NG)/(N-1)  approaches  (Pg.-1)G]^  in 
probability.  If  Pg.  = 1,  this  would  suggest  that  again  the 

distribution  of  Sj^  may  be  approximated  by  that  of  NGj.  From  the 

i 

! Central  Limit  Theorem, 

1/2 

Sj^  = NG^  + (N-l)(p^-l)Gj  + (N-1)  a^GiU^ 
where  Uj^  has  approximately  a standard  normal  distribution  for  large 
N.  Again,  the  goodness  of  the  approximation  depends  on  a^. 

Case  IV:  ~ ^j-1  3^1* 


This  looks  somewhat  like  Case  II.  In  fact,  it  is  quite  different. 
In  this  case: 


is  asymptotically  standard  normal.  One  could  again  say  that  Sj^  is 

approximated  by  when  = 0,  but  the  approximation  may  not  be 

1/2 

good.  In  cases  II  and  III,  was  multiplied  by  a factor  (N-1) 

3/2 

Here  , however,  is  multiplied  by  a factor  approximately  (N-1)  , 

which  is  larger  and  will  require  much  smaller  for  the  approximation 
to  be  plausible. 

There  is,  however,  another  feature  of  this  model  that  is  rather 
intriguing.  Namely,  it  can  bs  shown  that  (Appendix); 

Sn  - Nyg  - N(N-l)yg/2 
[N^Og^  + N(N-l)  (2N-l)ag,^/6]^/2 

has  an  asymptotic  standard  normal  distribution.  Note  that  appears 
nowhere  in  this  expression,  only  the  first  two  moments  of  Gj  appear 
in  the  expression  and  this  result  holds  whatever  the  distribution  of 
Gj,  provided  Og  is  finite.  Such  a result  is  very  nice  from  a 
statistical  viewpoint.  However,  we  shall  not  pursue  its  implications 
in  this  report. 

4.  THE  ROBUSTNESS  OF  RACICOT'S  MODEL 

As  previously  mentioned,  Racicot  assumed  that  the  Gj's  remain 
constant  throughout  the  life  of  the  tube.  However,  he  assumed  that 
the  material  and  design  parameters  entering  into  G varied  from  tube 
to  tube.  These  assumptions  will  be  adopted  here.  Equation  (6)  then 

^In  fact  similar  results  not  involving  Gj^  or  its  moments  in  any  way 
can  also  be  derived.  See  Appendix. 
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holds  for  any  particular  tube,  but  G,  k and  b^  will  vary  from  tube  to 
tube.  Equation  (5)  expresses  G as  a function  of  various  material 
parameters.  Each  of  these  parameters  is  assumed  to  vary. 

Fatigue  failure  occurs  when  the  crack  either  penetrates  through 
the  tube  or  reaches  a critical  depth  where  unstable  growth  occurs. 

The  critical  depth  b^  is  given  by: 

b = ^ (12) 

c TT  aS 

where  A is  an  empirical  constant  that  accounts  for  differences  in 
crack  shape  and  in  the  specimens  used  to  determine  Kj^.. 

Fatigue  failure  then  occurs  on  the  earliest  cycle  N such  that 
where  B is  the  wall  thickness  of  the  tube.  Since  N 
is  usually  large,  we  can  ignore  the  fact  that  N must  be  an  integer. 

This  gives 

N = (bo'^'-b"’')  (13) 

where  b = min(bg,B). 

If  the  probability  distributions  of  bg  and  of  each  of  the  model 
parameters  were  known,  then  it  would  be  relatively  simple  to  approximate 

the  distribution  of  N through  Monte  Carlo  simulation.  However,  none 
of  these  distributions  are  known.  Racicot  therefore  assumed  that  the 
distributions  of  each  of  the  relevant  random  variables  followed  some 
distribution  of  known  mathematical  form  and  then  estimated  the 
parameters  of  the  distributions  from  test  data  gathered  by  Davidson 
and  others.  He  then  generated  the  distribution  of  N by  Monte  Carlo 
simulation  using  equation  (13) . 
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Certain  methods  used  by  Racicot  are,  however,  open  to  criticism 


and  could  possibly  weaken  his  conclusions. 

First,  Racicot  assumed  that  each  of  the  model  parameters  had  the 

s^lme  type  of  distribution.  Specifically,  he  considered  3 cases: 

all  of  the  parameters  had  normal  distributions;  all  had  Weibull 

distributions;  and  all  had  lognormal  distributions.  It  may  in  fact 

be  the  case  that  all  the  model  parameters  have  the  same  type  of 

distributions,  at  least  approximately.  But  there  is  at  this  time  no 

real  reason  to  believe  so.  Even  if  they  did  have  the  same  type  of 

distribution,  there  is  no  reason  to  assume  that  this  type  would  be 

either  normal,  lognormal  or  Weibull. 

Second,  Racicot  compared  his  Monte-Carlo  distribution  with  various 

possible  candidate  distributions  by  using  the  Kolmogorov- Smirnov  (KS) 

5 . 

statistic  with  the  parameters  of  the  candidate  distribution  estimated 
^ 4 

by  the  moment -matching  method  . However,  for  the  lognormal  and 
Weibull  distributions,  he  moment -matching  method  is  not  the  best 
method  to  use.  Further,  the  KS  statistic  is  known  to  be  inferior  in 
some  ways  to  other  goodness-of-fit  statistics^.  How  would  Racicot 's 
conclusions  fare  under  better  methods  of  estimation  and  other  goodness- 
of-fit  tests? 

'^Cramer,  H.,  Mathematical  Methods  of  Statistics,  Princeton  University 
Press,  Princeton,  NJ,  1946,  p.  497. 

^Lindgren,  B.  W.,  Statistical  Theory,  The  Macmillan  Company,  New  York, 
NY,  1970,  p.  329. 

^Stephens,  M.  A.,  "EOF  Statistics  for  Goodness  of  Fit  and  Some 
Comparisons,"  Journal  of  the  American  Statistical  Association,  69, 
(September  1974),  pp.  730-737. 


Third,  other  candidate  distributions  could  be  considered. 

Fourth,  Racicot  assumed  that  most  of  the  model  parameters  were 
statistically  independent.  Again,  there  is  no  real  reason  to  assume 
that  this  is  so. 

The  next  sections  of  this  report  will  describe  the  efforts  of  the 
author  to  fill  in  these  gaps. 

5.  INDEPENDENT  MODEL  PARAMETERS  WITH  RANDOMLY  SELECTED  DISTRIBUTIONS 
Racicot  used  data  from  nine  105mm  M137A1  tubes  tested  by  Throop 
and  others  as  a base  for  estimating  the  distributions  of  the  model 
parameters.  The  present  author  will  continue  to  use  these  data. 

«,  E,  bjj  and  m are  assumed  to  be  random  variables  with  means 

and  standard  deviations  as  given  in  Table  1. 

TABLE  1.  MEANS  AND  STANDARD  DEVIATIONS  OF  MODEL  PARAMETERS  FOR 


105MM  M137A1 

TUBES 

Parameter 

Mean 

Standard  Deviation 

a 

0.877 

0.050 

Kjj,  (KSl/in) 

103 

11.5 

E (KSI) 

30,000 

300 

bo  (in) 

0.02 

0.001 

m 

3.5 

0.1 
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Oy.  was  given  by  the  empirical  relationship: 

Oy  = 334  - 1.39  Kjj,.  The  other  parameters  (S,A,B,C)  were  assumed 
constant,  a,  Kj^.,  E,  Eq  and  m were  assumed  to  be  statistically 
independent . 

The  distribution  of  N,  the  fatigue  failure  time,  was  then 
approximated  by  computer  simulation.  For  each  run  of  the  simulation 
program  through  the  computer,  the  program  itself  chose,  at  random, 
a type  of  distribution  of  each  of  the  five  random  model  parameters. 

The  choice  of  a distribution  type  for  each  of  the  model  parameter  was 
independent  of  the  choice  of  type  for  any  of  the  others.  The  ration- 
ale behind  this  method  is:  The  author  wanted  to  check  that  Racicot's 
assumption  of  the  same  distribution  type  for  each  of  the  five  parameters 
did  not  somehow  bias  his  results.  However,  it  is  clearly  impossible 
to  try  all  possible  different  combinations  of  types.  By  instructing 
the  computer  to  select  distribution  types  at  random,  the  author  hoped 
to  obtain  a representative  sampling  of  possible  combinations.  The 
possible  distribution  types  available  for  selection  were:  normal, 
lognormal,  uniform,  extreme  value,  Weibull  and  2-parameter  exponential. 

Once  the  distribution  type  was  selected,  the  actual  parameters 
of  the  distribution  were  estimated  by  the  method  of  moments  from  the 
means  and  standard  deviations  given  above. 

Ten  thousand  replicates  of  N were  then  generated  for  each  run  by 
generating  ten  thousand  sets  of  replicates  of  o,  Kj^.,  E,  b^  and  m 
and  calculating  a corresponding  N for  each  set  of  replicates  using 
equation  (13).  Because  of  computer  space  limitations,  it  was 
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necessary  to  round  each  N to  the  nearest  10  cycles;  however,  this 
should  have  no  significant  effect. 


The  distribution  of  N generated  by  this  method  was  then  compared 
to  each  of  nine  candidate  parametric  distribution  families.  For 
uniformity,  we  will  characterize  each  of  the  statistical  parameters 
as  either  a location  parameter  denoted  by  n.  a shape  parameter 
denoted  by  $ or  a scale  parameter  denoted  by  0.  Also,  ♦(x)  represents 
the  standard  normal  cumulative  distribution  function,  and  F(x) 
represents  the  cumulative  distribution  function  (cdf) . The  nine 
families  are: 


Normal 


Lognormal 


F(x)  = $(21111) 

0 

FCx)  = Inf)] 


lixtreme  value 

F(x)  = l-exp[-e  ] 


Weibull 


F(x)  = l-exp[-(x/0)  ] 


1 - parameter  Exponential 

F(x)  = 1-exp (-x/0) 

2- parameter  Exponential 

F(x)  = 1-exp  (13^) 

Double-Exponential  (Fisher-Tippett  Type  I) 
F(x)  = exp[-e^^"^^^®] 


Inverse-Weibull 


Bimbaum-Saunders 


F(x)  = exp[-(0/x)n 

F(x,  = 


all  x 

X > 0 

all  X 

X > 0 

X > 0 

X > ri 

all  X 

X > 0 

X > 0 
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For  all  of  these,  F(x)  = 0,  except  where  otherwise  noted. 

The  3-par2imeter  versions  of  the  Weibull,  lognormal  and  inverse- 
Weibull  will  not  be  considered  here.  These  distributions  are  very 
difficult  to  deal  with  in  practice  with  the  limited  data  found  in 
fatigue-life  problems.  The  whole  point  of  this  exercise  is  to  give 
the  applied  statistician  some  method  he  can  use  in  real-life  problems. 
At  the  present  time,  the  3-parameter  versions  simply  do  not  provide 
such  a method. 

It  is  also  generally  true  that  the  more  parameters  one  adds  to  a 
distribution,  the  better  one  can  fit  arbitrary  sets  of  data  to  it, 
even  if  the  distribution  has  little  justification  on  empirical 
grounds.  Multi-parameter  Weibull  distributions  have  given  good  fits 
to  some  fatigue  data  in  the  past  probably  because  they  had  enough 
parameters  to  give  a wide  variety  of  shapes  By  limiting  ourselves 
to  two-parameter  distributions,  we  will  hopefully  arrive  at  something 
that  has  some  relation  to  underlying  physical  process  and  avoid 
labeling  a distribution  as  good  fitting  only  because  it  has  enough 
parameters  to  smooth  out  the  rough  spots. 

For  each  of  the  families,  approximation  of  the  statistical 
parameters  (ti,0,  and  B)  used  the  moment -matching  method  for  the 
normal,  extreme-value,  1-parameter  exponential,  double-exponential 
and  Birnbaum-Saunders.  Moment -matching  on  the  logarithms  of  the 

^Little,  R.  E.,  and  Jebe,  E.  H.,  The  Statistical  Design  of  Fatigue 
Experiments,  New  York,  John  Wiley  S Sons,  1975,  p.  221. 
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fatigue  lives  was  used  for  estimation  for  the  lognormal,  Weibull  and 
inverse- Weibull.  The  2-parameter  exponential  used  minimum-variance 
unbiased  estimators. 

When  estimators  of  the  statistical  parameters  are  substituted  for 
the  actual  parameters  in  the  appropriate  expression  for  F{x),  one 
obtains  an  estimate  F(x)  of  F(x).  The  function  F(x)  can  then  be 
compared  to  the  Monte  Carlo  distribution  of  N to  test  goodness  of  fit. 
We  have  used  three  different  goodness-of-fit  statistics.  They  are  as 
follows: 


Kolmogorov-Smirnov  (KS) 


D = max  [ i - Z^,  ] 

l<i<M  M 1 1 M 


Cramer-von-Mises  (CVM) 

W^  = — i-  + 

IZM^ 

Anderson- Darling  (AD) 


- ^ (Zi  - 
M i=l  ^ 


2i-U2 
2M  ^ 


M Zi 

I (2i-l)ln(-  ■■  >— -) 
i=l  I-Zm+1-1 


where  M = 10,000 


= the  ith  smallest  fatigue  life  from  the  Monte  Carlo 
distribution  (1  £ i £ 10,000) 

Z.  = F(Ni) 

For  all  three  of  these  goodness-of-fit  statistics,  the  smaller 
the  value,  the  better  the  fit. 

This  whole  procedure  was  repeated  seven  times,  giving  seven 
different  possible  distributions  for  the  fatigue  lives.  Details  are 
given  here  for  three  of  them,  which  are  representative  of  the  type  of 
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results  obtained.  Table  2 shows  the  distribution  types  selected  for 
each  of  the  five  varying  model  parameters.  Tables  3,  4 and  5 show 
the  results  of  the  goodness-of-fit  tests. 

For  run  #1,  lognormal  gives  the  best  fit,  followed  very  closely 
by  Birnbaum-Saunders . These  results  are  reversed  for  run  #2.  However, 
it  should  be  noted  that  the  cdf's  of  the  lognormal  and  the  Birnbaum- 
Saunders  are  so  close  as  to  be  practically  indistinguishable.  However, 
the  fit  is  not  of  the  best.  In  fact,  all  three  goodness-of-fit 
statistics  for  the  lognormal  are  significantly  too  large  at  a 
significance  level  of  15%  or  smaller 

For  run  #3,  the  inverse-Weibull  gives  an  excellent  fit. 
Unfortunately,  tables  for  testing  the  goodness-of-fit  statistics  in 
this  case  are  not  available. 

Figures  1,  2 and  3 show  histograms  of  the  Monte  Carlo  fatigue 
life  simulations.  Figures  4,  5 and  6 show  the  cdf  of  the  Monte  Carlo 
distribution  along  with  the  lognormal,  Weibull,  inverse-Weibull  and 
2-parameter  exponential  distributions  used  to  approximate  it. 

6.  CORRELATED  MODEL  PARAMETERS 

Up  to  now,  the  five  model  parameters  a,  Kic,  E,  bg  and  m have  been 
considered  statistically  independent.  There  is  no  real  reason  for 
believing  that  they  actually  are. 


^Stephens,  M.  A.,  "EDF  Statistics  for  Goodness  of  Fit  and  Some  Compar- 
isons," Journal  of  the  American  Statistical  Association,  69,  (Septem- 
ber 1974),  pp.  730-737.  (Stephen's  tables  for  Case  2 normal  apply 
unchanged  to  the  lognormal  distribution  with  the  type  of  estimators 
we  have  been  using) . 
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TABLE  2.  DISTRIBUTION  TYPES  OF  THE  MODEL  PARAMETERS 


1 


Run  #I 

Run  #2 

Run  #3 

a 

uniform 

uniform 

extreme-value 

•^Ic 

normal 

normal 

uniform 

E 

uniform 

lognormal 

lognormal 

bo 

uniform 

uniform 

extreme-value 

m 

extreme-value 

normal 

lognormal 

I 
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TABLE  3.  GOODNESS-OF-FIT  RESULTS 
Run  #1 


Distribution 


Parameters 


Goodness-of-Fit  Statistics 


Type 

Location 

Shape 

Scale 

KS 

CYMCxlO**! 

AD 

Normal 

11737 

- 

2570 

0.072 

16.6 

10.6 

Lognormal 

- 

4.594 

11463 

0.051 

9.5 

6.2 

Extreme-Value 

12894 

- 

2004 

0.13 

61.9 

48.9 

Weibull 

- 

5.893 

12642 

0.10 

38.7 

29.5 

Exponential 

1-parameter 

- 

- 

11738 

0.47 

594.7 

282.5 

Exponential 

2-parameter 

5630 

- 

6108 

0.30 

287.3 

149.9 

Double- 

Exponential 

10581 

- 

2004 

0.063 

16.5 

10.8 

Inverse- 

Weibull 

- 

5.893 

10393 

0.087 

28.0 

20.2 

Bimbaiim- 

Saunders 

0.2178 

11466 

0.052 

9.6 

6.2 

TABLE  4.  GOODNESS-OF-FIT  RESULTS 


Run  #2 

Distribution  Parameters  Goodness-of-Fit  Statistics 


Type Location 


Normal 

11689 

Lognormal 

- 

Extreme-Value 

12960 

Weibull 

- 

Exponential 

1 -parameter 

- 

Exponential 

2-parameter 

5310 

Double- 

Exponential 

10418 

Inverse- 

Weibull 

- 

Bimbaum- 

Saunders 

Shape 

Scale 

KS 

- 

2825 

0.062 

4.182 

11359 

0.035 

- 

2202 

0.129 

5.363 

12650 

0.090 

- 

11689 

0.445 

- 

6379 

0.281 

- 

2202 

0.044 

5.363 

10200 

0.080 

0.2401 

11361 

0.035 

CVMCxlO  1 AP 


12.5 

7.9 

3.8 

2.4 

58.4 

>10® 

29.9 

23.9 

562.6 

268.5 

268.8 

140.2 

8.1 

5.22 

21.9 

16.69 

3.7 

2.3 
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« TABLE  5.  GOODNESS-OF-FIT  RESULTS 


Run 

#3 

Distribution 

Parameters 

Goodness 

-of-Fit 

Statistics 

! Type 

Location  Shape 

Scale 

KS 

CVMCxlO** 

) AD 

1 : 

1 Normal 

12627 

5834 

0.153 

99.1 

>10® 

Lognormal 

2.854 

11768 

0.068 

17.0 

10.5 

Extreme-Value 

15253 

4549 

0.220 

179.1 

>10® 

' Weibull 

! 

3.660 

8324 

0.137 

68.9 

>10® 

; Double 

Exponential 

10001 

4549 

0.120 

44.6 

30.1 

1 Exponential 

1 -parameter 

- 

12627 

0.399 

425.5 

206.1 

Exponential 

2-parameter 

4950 

7677 

0.229 

156.4 

85.9 

i Inverse- 

i Weibull 

3.660 

10051 

0.011 

0.2 

0.2 

Bimbaum- 

Saunders 

0.4545 

11445 

0.117 

40.9 

28.2 

! 


20 


t . . . 

_ . 



1 , . 

} 1 

, , - - 

...  - 

L.  . 

' • “ I 

r 

: I ; : 

I _ 

' I 

J 

o 

! ■ 

I 

; 

. . 

^___ 

H 

- J 

I 

J 

I 

L 

I 

' 

Tl 

• - 

’ 1 

J 

. 1 

Two  more  runs  were  performed  as  described  in  the  last  section. 

However,  in  these  two  runs,  all  of  the  five  model  parameters  were  ; 

taken  to  be  normal.  On  each  run,  the  program  selected  at  random  a 
correlation  coefficient  for  each  pair  of  the  five  model  parameters. 

1 

The  results  for  these  two  runs  (Labeled  NCI  and  NC2)  are  shown  in  j 

Tables  6 through  8 and  in  Figures  7 through  10.  For  NCI,  the  lognormal  j 

gives  an  excellent  fit.  None  of  the  three  goodness-of-flt  statistics  | 

are  significant  at  any  level  below  15%.  For  NC2,  the  double 
exponential  gives  the  best  fit,  followed  along  by  the  lognormal.  The 
lognormal  does  not  give  too  bad  a fit.  The  KS  statistic  is  significant 
at  15%  and  below,  but  neither  the  CVM  nor  the  AD  statistics  are 
significant. 

7.  THE  EFFECTS  OF  THE  MODEL  PARAMETER  a 

Racicot  concluded  that  the  model  parameter  a contributed  more  than 
any  of  the  other  parameters  to  the  variability  in  fatigue  life.  A 
stronger  statement  can  be  made,  a is  not  only  the  major  source  of 
variability  in  fatigue  life,  but  it  also  is  the  major  determinant  of 
the  very  shape  of  the  fatigue  life  distribution.  It  would  be  an 
exaggeration  to  say  that  E,  Kj^.,  bQ  and  m are  all  second-order  effects 
compared  to  a;  however,  it  would  not  be  a gross  exaggeration. 

Figures  11,  12  and  13  show  the  cdf  of  fatigue  life  obtained  when 
only  a is  assumed  to  vary  from  tube  to  tube.  (The  other  model 
parameters  are  kept  constant  at  their  mean  values.)  The  dotted  lines 
show  the  cdf's  obtained  when  all  five  of  the  model  parameters  a,  Kj^., 

E,  bg  and  m are  assumed  to  vary  from  tube  to  tube.  Figure  11  shows 


TABLE  6.  CORRELATION  COEFFICIENTS  FOR  THE  MODEL  PARAMETERS 


lei  Parameters 

NCI 

NC2 

ot,  Kjg 

0.719 

-0.406 

a,  E 

-0.498 

0.639 

a,  bo 

0.111 

0.408 

a,  m 

0.097 

-0.610 

Kic,  E 

-0.937 

0.391 

bo 

0.175 

0.437 

Kjg,  m 

0.339 

-0.415 

E,  bo 

-0.357 

0.609 

E,  m 

-0.372 

-0.999 

bo.  ™ 

0.462 

-0.584 

TABLE  7.  GOODNESS-OF-FIT  RESULTS 


Run  NCI 


Distribution 

Type 

Parameters 
Location  Shape 

Scale 

Goodness-of-Fit 
KS  CVMfxlO'*') 

Statistics 

AD 

Normal 

11620 

- 

2262 

0.043 

6.9 

4.5 

Lognormal 

- 

5.224 

11408 

0.007 

0.08 

0.08 

Extreme-Value 

12638 

- 

1763 

0.114 

47.5 

>10® 

Weibull 

- 

6.700 

12435 

0.075 

21.0 

>10® 

Exponential 

1- parameter 

- 

- 

11620 

0.469 

649.0 

305.2 

Exponential 

2-parameter 

5169 

- 

6451 

0.349 

406.0 

200.5 

Double 

Exponential 

10602 

- 

1763 

0.033 

3.52 

2.7 

Inverse- 

Weibull 

- 

6.700 

10466 

0.071 

18.1 

14.9 

Birnbaum- 

Saunders 

0.1938 

11406 

0.009 

0.2 

0.2 
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TABLE  8.  GOODNESS-OF-FIT  RESULTS 
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Run 

NC2 

Distribution 

Parameters 

Goodness-of-Fit 

Statist: 

Type 

Location 

Shape 

Scale 

KS 

CVMfxlO**! 

AD 

Normal 

11879 

- 

3210 

0.068 

15.9 

>10® 

Lognormal 

- 

3.867 

11482 

0.016 

0.5 

0.3 

Extreme-Value 

13323 

- 

2502 

0.138 

66.9 

>10* 

Weibull 

- 

4.960 

12899 

0.084 

25.6 

>10* 

Exponential 

1 -parameter 

- 

- 

11879 

0.426 

542.0 

258.8 

Exponential 

2 -parameter 

4389 

- 

7489 

0.310 

322.0 

162.4 

Double- 

Exponential 

10434 

- 

2502 

0.010 

0.3 

0.2 

Inverse- 

Weibull 

- 

4.960 

10220 

0.061 

14.7 

12.7 

Bimbainn- 

Saunders 

0.2681 

11466 

0.018 

0.8 

0.7 

Frequency  histogram  - Run  NC2 


the  cases  when  a is  assumed  to  be  uniform;  Figure  12,  when  a is 

normal;  and  Figure  13,  when  a is  extreme-value.  Except  for  one  of  the 

dotted  curves  on  Figure  13  (corresponding,  incidentally,  to  run  #3 

discussed  above),  the  dotted  curves  have  approximately  the  same  shape 

as  the  solid  curves  and  tend  to  remain  fairly  close  together. 

8.  SOME  COMMENTS  ON  THE  PROSCHAN-SETHURAMAN  MODEL 

2 

Proschan  and  Sethuraman  have  developed  a model  for  gun  tube 
fatigue  also  based  on  the  Throop  equation.  They  assume,  as  did 
Racicot,  that  none  of  the  material  parameters  vary  from  cycle  to 
cycle.  They  also  assume,  however,  that  none  of  the  material  properties 
vary  from  tube  to  tube.  In  the  notation  of  this  report,  they  assume 
that  G is  constant  from  tube  to  tube.  They  do  assume  that  bp  varies 
from  tube  to  tube. 

Further,  Proschan  and  Sethuraman  assume  that 

a.  A large  nimber  of  initial  cracks  are  present  and  the  largest 
one  present  ultimately  causes  the  failure. 

b.  The  number  of  initial  cracks  has  a Poisson  distribution  and 
the  depths  of  the  initial  cracks  are  independent  and  identically 
distributed. 

c.  There  exists  some  constant  y such  that  all  initial  crack 
lengths  are  no  larger  than  y,  with  probability  1. 

2Proschan,  F.,  and  Sethuraman,  J.,  "A  Probabilistic  Model  for  Initial 
Crack  Size  and  Fatigue  Life  of  Gun  Barrels,"  Florida  State  University 
Technical  Report  No.  M394,  The  Florida  State  University,  Tallahassee, 
Florida,  December  1976. 
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The  present  author  finds  assumptions  a,  b and  c reasonable  and 
cannot  quibble  with  them.  He  can  (and  will)  quibble  with  the 
assumption  that  all  material  properties  are  constant  from  tube  to 
tube.  From  the  results  of  the  last  section,  it  is  clear  that  the 
variation  in  the  material  property  a Ccinnot  be  ignored. 

Proschan  and  Sethuraman  state  without  proof  that  "it  can  be  shown 
from  the  theory  of  limiting  distributions  for  the  maximum  (or  minimum) 
of  large  numbers  of  random  variables  that  the  distribution  of 

2 

[fatigue  life]  can  be  approximated  by  a shifted  Weibull  distribution"  • 

It  may  be  interesting  to  note  that  under  their  assumptions  the  fatigue  life 
will  approach  a constant  in  probability  (and  hence  also  in  distribution) 
as  the  average  number  of  initial  cracks  becomes  large.  This  is  easily 
seen  as  follows:  Let  M be  the  number  of  initial  cracks  and  let 
Dj,..  be  their  depths.  Let  D = max  (Dj^  ,•  • • » D]q)  . By  assumption  a, 
fatigue  life  is  a continuous  function  of  D.  So  it  suffices  to  show 
that  D approaches  a constant  in  probability.  By  assumption  b, 
have  a common  distribution.  Denote  its  cdf  by  F.  By 
assumption  C,  there  exists  some  y such  that  F(x)  = 1 if  x > y and 
F(x)  < 1 if  X < y. 

Now,  for  any  x and  M ^ 1, 

Pr(D  < x|M)  = Pr(Dj  1 x,-  • -.D^  _<  x|M)  = F^(x) 

2 

Proschan,  F.,  and  Sethuraman,  J.,  "A  Probabilistic  Model  for  Initial 
Crack  Size  and  Fatigue  Life  of  Gun  Barrels,"  Florida  State  University 
Technical  Report  No.  M394,  The  Florida  State  University,  Tallahassee, 
Florida,  December  1976,  p.  6. 


r 


M has  a Poisson  distribution  with  mean  X.  (X  is  then  the  average 
number  of  initial  cracks  in  a tube.)  Then 

^ Pr(D  < xIm  >1) 

Pr(D  _<  x)  = — 1-^= — ! — 

Pr(M  > 1) 

00 

j;e"^X^F^(x)/M  ! 

= M^l 

I e'H^M'. 

M=1 

^ e>^FCx).i 


.-X[l-F(x)]  -X 


Now  let  e > 0.  Then  Pr(D  ^ y G)  = 0.  So 
Pr(|D-Yl  ^ e)  = Pr(D  £ y-e) 

e-X[l-F(Y-e)].e-^ 


But,  Y-e  < Y and  so  l-F(Y-e)  > 0.  Then  as  X -►  »,  it  follows  from 
equation  (14)  that  Pr(|D-Y|  ^ e)  -*■  0. 

Hence,  D -►  y in  probability  as  X -*■  «>. 

9.  CONCLUSIONS 

Racicot  concluded  from  his  research  that  the  lognormal  gave  the 
best  fit  to  the  fatigue  life  distribution.  The  present  study  shows 
that  Racicot 's  conclusion  holds  up  fairly  well  under  conditions 
broader  than  he  considered.  The  lognormal  (and  its  close  ally  the 
Bimbaum-Saunders)  almost  always  gave  the  best  fit  or  ran  a close 
second.  However,  it  must  be  stressed  that  the  lognormal  can  only  be 
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an  approximation.  This  is  no  reason  to  assume  that  the  fatigue  life 
distribution  ^ lognormal.  In  fact,  assumptions  of  lognormality 
could  be  rejected  in  most  cases. 

What  is  significant  is  that  the  lognormal  performed  much  better 
than  its  principal  rival,  the  Weibull,  which  has  often  been  used  in 
fatigue  life  studies.  The  two-parameter  exponential,  suggested  by 
Proschan  and  Sethuraman,  is  not  even  close.  In  this  regard,  it  is 
interesting  to  note  that  the  distributions  such  as  the  inverse-Weibull, 
which  to  the  author's  knowledge  has  never  been  mentioned  as  a 
possibility  for  approximating  the  fatigue  life  distribution,  did 

better  than  the  much  mentioned  Weibull.  i 

I 

The  author  would  have  been  quite  delighted  to  be  able  to  announce  l 

i 

that  a lognormal  (or  Weibull,  or  inverse-Weibull,  etc.)  was  so  close  to  i 

the  fatigue  life  distribution  that  it  could  be  used  with  no  further  j 

qualms.  It  does  not  appear  at  this  time  that  any  such  announcement  can 
be  made.  There  would  also  seem  to  be  little  use  quibbling  as  to 
whether  the  fatigue  life  distribution  actually  is  lognormal  or  Weibull 
or  exponential.  It  is  most  likely  none  of  these. 

At  this  point,  a more  useful  approach  would  be  to  drop  the  attempt 
to  find  out  what  the  fatigue  life  distribution  really  is,  and  instead 
to  discover  if  the  procedures  derived  from  the  candidate  distributions 
give  good  results  even  though  the  actual  fatigue  life  distribution  is 
none  of  the  candidates.  In  other  words,  we  should  examine  such  matters 
as  to  whether  a confidenced  safe  life  computed  under  an  assumption  of 
(say)  lognormality  for  the  fatigue  lives  gives  a reasonable  approximation 

i 

i 
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to  a true  confidenced  safe  life.  Answering  this  type  of  question  could 
ultimately  be  more  fruitful  than  engaging  over  arguments  over  the  true 
nature  of  the  fatigue  life  distribution. 
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APPENDIX 

MATHEMATICAL  DETAILS 

This  Appendix  will  deal  with  the  mathematical  details  of  the 
assertions  made  in  Section  3,  Case  IV.  For  this  we  will  need  a slight 
generalization  of  the  Central  Limit  Theorem. 

Theorem;  Let  (Xjj)  be  a sequence  of  independent,  identically 
distributed  random  variables  with  Ex,^  = y and  Var(Xj^)  = where 
0 < a < “.  For  each  n,  let  “nr‘’’®4in  constants.  Let  (x^)  and 
(Pu)  be  sequences  of  constants  such  that; 
n 


a)  ( I Onjl^-^n) 

Tn  j=i 


b) 


max  , “nj  , 

l<j<n  I I 0 

n 


as  n 00 


as  n -+■  00 


as  n “ 


1 " 

Then  ( ^ bas  an  asymptotic  standard  normal 

cf^n  j = i ^ 

distribution. 

Proof;  This  theorem  is  proved  essentially  the  same  way  as  the 

4 

standard  Central  Limit  Theorem  , which  is  the  special  case  = y. 
Tjj  = l/i/n,  Ujjj  = 1/n. 

Let  <()(u)  be  the  common  characteristic  function  of  the  random 

variables  (x^j  - y)/o.  Let  be  the  characteristic  function  of 

1 ” 


^Cramer,  H. , Mathematical  Methods  of  Statistics,  Princeton  University 
Press,  Princeton,  NJ,  1946,  p.  213.  (The  proof  of  the  theorem  is 
practically  the  same  as  Cramer's  proof  of  the  Central  Limit  Theorem, 
so  many  details  have  been  omitted.) 
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Then 


Now  assume  the  set-up  as  in  Case  IV,  Section  3.  There  we  had 
N-1 

Sn  - NGj  - ^ (N-j)ej  where  the  ej's  are  independent  and  identically 

2 

distributed  with  mean  vij.  and  variance  . The  theorem  can  now  be 

applied  with  aj,j  = n + 1 - j for  1 £ j ^ n» 

_ _ The  formulas® 

^ 2 


n(n-t-l)  (2n-*-l) 
6 


j=l  ^ j=l 

may  be  used  to  verify  that  the  conditions  of  the  theorem  hold.  The 

conclusion  is  that 

Sjj  - NGj  - N(N-l)y6/2 

[N(N-n  (2N-2)aeV6]^^^ 

has  an  asymptotic  standard  normal  distribution. 

Now, 

N-1 

E(Sj^)  = Nyg  + ^ (N-j)yg. 

j=l 


NUg 


(16) 


If  G,  is  independent  of  the  e's,  then 

N-1 

Var(S^)  = N^og*  + [ (N-j)^a^^ 

j = l 

= . N(N-1).(2.N-1J.  2 

G ^ e 


(17) 


N^a,  ^ 


Let  T 


G_  ^ N(N-1)(2N-1)  . Then  N^^^/tn  ^ /3 


N 


®Selby,  S.,  Standard  Mathematical  Tables  C14th  Edition),  The  Chemical 
Rubber  Co.,  Cleveland,  OH,  1964,  p.  390. 
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r 


Then 


0. 


N . . . N _ 1 

- lif;)  -*■  0 in  probability,  since  — = , 

a^TN  G ^\/2  Tn 

Application  of  the  theorem  gives  that  — - — C I (N-j)e^  - y ) 

OeTn  j=i  J 2 ^ 

has  an  asymptotic  standard  normal  distribution.  Then 

Sn  - Npg  - N(N-l)yg 


N-1 


= -N-  (Gi-yg)  +-i_  ( I (N-j)e. 

a^TN  j=l 

g 

also  has  an  asymptotic  standard  normal  distribution  . This  last 
result  does  not  require  independence  of  from  the  Ej's,  although 
equation  (17)  would  not  necessarily  hold  without  such  independence. 

It  is  also  possible  to  obtain  asymptotic  results  not  involving  Gj 
or  its  moments  in  any  way.  In  fact  if  (tj^)  is  any  sequence  such  that 
/3,  then 


a) 


NGl 


0 in  probability; 


N-1 


b)  — i \ (N-j)e.  - y ) has  an  asymptotic  standard  normal 

Vn  ^ 2 ^ 

distribution  (by  the  Theorem) . 

Then 

ag.Tfj  Oj-Tj^  Og-Tf^  j = l 2 

has  an  asymptotic  standard  normal  distribution. 


'Rao,  C.  R. , Linear  Statistical  Inference  and  Its  Applications, 
John  Wiley  5 Sons,  Inc.  New  York,  NY,  1965,  p.  102,  Theorem  (x) . 
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